J Arid Land 
https://doi.org/10.1007/s40333-018-0006-1 


Characteristics of soil organic carbon and total 
nitrogen under various grassland types along a 


transect in a mountain-basin system in Xinjiang, 
China 


BI Xu', LI Bo'’", NAN Bo', FAN Yao!, FU Qi!, ZHANG Xinshi!? 


' Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China; 
? Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China 


Abstract: Soil organic carbon (SOC) and soil total nitrogen (STN) in arid regions are important 
components of global C and the N cycles, and their response to climate change will have important 
implications for both ecosystem processes and global climate feedbacks. Grassland ecosystems of Funyun 
County in the southern foot of the Altay Mountains are characterized by complex topography, suggesting 
large variability in the spatial distribution of SOC and STN. However, there has been little investigation of 
SOC and STN on grasslands in arid regions with a mountain-basin structure. Therefore, we investigated 
the characteristics of SOC and STN in different grassland types in a mountain-basin system at the 
southern foot of the Altai Mountains, north of the Junggar Basin in China, and explored their potential 
influencing factors and relationships with meteorological factors and soil properties. We found that the 
concentrations and storages of SOC and STN varied significantly with grassland type, and showed a 
decreasing trend along a decreasing elevation gradient in alpine meadow, mountain meadow, temperate 
typical steppe, temperate steppe desert, and temperate steppe desert. In addition, the SOC and STN 
concentrations decreased with depth, except in the temperate desert steppe. According to Pearson's 
correlation values and redundancy analysis, the mean annual precipitation, soil moisture content and soil 
available N concentration were significantly positively correlated with the SOC and STN concentrations. 
In contrast, the mean annual temperature, pH, and soil bulk density were significantly and negatively 
correlated with the SOC and STN concentrations. The mean annual precipitation and mean annual 
temperature were the primary factors related to the SOC and STN concentrations. The distributions of 
the SOC and STN concentrations were highly regulated by the elevation-induced differences in 
meteorological factors. Mean annual precipitation and mean annual temperature together explained 
97.85% and 98.38% of the overall variations in the SOC and STN concentrations, respectively, at soil 
depth of 0—40 cm, with precipitation making the greatest contribution. Our results provide a basis for 
estimating and predicting SOC and STN concentrations in grasslands in arid regions with a 
mountain-basin structure. 
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1 Introduction 


The carbon (C) and nitrogen (N) cycles of terrestrial ecosystems have received increasing 
attention worldwide in recent decades due to their sizable roles in global climate change and 
ecosystem sustainability (Fu et al., 2010). Soils represent the largest terrestrial C pool and 
contribute to nitrous oxide (N2O) emissions (Saikawa et al., 2014; Nocita et al., 2015). Soil N can 
also deeply impact ecosystem productivity and C sequestration (Dixon et al., 1994; Reich et al., 
2006). Therefore, a better understanding of the characteristics and impact factors underlying soil 
C and N at both regional and global scales is critical for evaluating C and N cycles and for 
predicting the feedback of soil C and N to global climate change (Jobbágy and Jackson, 2000; 
Yang et al., 2007). As one of the major terrestrial ecosystems, grasslands in China cover 
approximately 40% of the total terrestrial surface and represent 13% of the total grasslands in the 
world. Grasslands are the most important terrestrial C pools in China, with 41.03 Pg C stored in 
the soil (Ni, 2002; Chen et al., 2016). During the past several decades, a number of studies on soil 
C and N, have focused on the grasslands at different scales in China, including the site scale (Yu 
et al., 2013; Wang et al., 2014; Wang et al., 2017), the transect scale (Zhou et al., 2002; Yang et 
al., 2009; Yang et al., 2010), and the country scale (Ni, 2002; Piao et al., 2005; Xie et al., 2007; 
Fang et al., 2010). However, most studies have assessed C and N storage values for specific 
grassland types, especially temperate grasslands (Conant and Paustian, 2002; Piao et al., 2007; 
Yang et al., 2010) and alpine grasslands (Li et al., 2014; Yan et al., 2015; Chen et al., 2017). 
Additionally, numerous studies have studied soil C and N in the grasslands of Inner Mongolia 
(Bai et al., 2004; He et al., 2014; Li et al., 2017) and the Qinghai-Tibetan Plateau (QTP) (Yang et 
al., 2008; Chang et al., 2014) in China. However, few studies have investigated soil organic 
carbon (SOC) and soil total nitrogen (STN) storages variations among different grassland types in 
the southern foot of the Altai Mountains due to limited soil surveys and the high spatial variability 
of soils in the arid region. 

Arid regions cover approximately 12.0% of the land area in the Earth (Lal, 2004). A 
distribution of alternating mountains and basins is the basic characteristic of the natural 
geography of the ecosystems in arid regions in Central Asia (Fu et al., 2017). Zhang (2001) 
defined this landform as a mountain-basin system (MBS) (Fig. 1). Therefore, the grasslands in 
arid regions with a MBS structure are an important component in the global C and N cycles, and 
their response to climate change will have important implications for both ecosystem processes 
and global climate feedbacks (Fang et al., 1996). 

Fuyun County lies in the arid region of Xinjiang, China, and it is a representative MBS, 
consisting of a vertical mountain vegetation gradient and a concentric annular vegetation gradient 
in a desert basin. Fuyun County represents an ideal region for studying soil C and N cycles in 
different grassland types and the corresponding feedback interactions with global climate change. 
Studying the SOC and STN in the grassland ecosystems of Fuyun County can improve our 
knowledge of grassland C and N cycles and provide guidance for assessing and predicting C and 
N storage in grasslands in arid region. Thus, we studied the characteristics of SOC and STN and 
their relationships with meteorological factors and soil properties in different grassland types in 
Fuyun County through field investigations and sampling. 

Soil C and N are often tightly coupled and are significantly impacted by many factors, such as 
climate (Dai and Huang, 2006; Mufioz-Rojas et al., 2017), elevation (Lemenih and Itanna, 2004; 
Tsui et al., 2013), vegetation (Fu et al., 2010; Roman-Sanchez et al., 2016) and soil properties 
(Sakin, 2012; Six and Paustian, 2014). Climate is commonly considered the most important factor 
regulating SOC and STN storages (Callesen et al., 2003;Wynn et al., 2006), because it affects the 
balance of C and N inputs from soil organic matter (SOM) and outputs through decomposition 
and mineralization in soils (Garten et al., 1999; Hevia et al., 2003). Precipitation affects SOC 
storage by constraining plant production and decomposition, especially in arid ecosystems (Wynn 
et al., 2006; Yang et al., 2007). Some studies have shown that higher temperatures tend to 
accelerate the microbial decomposition of soil organic matter, thereby causing C loss (Giardina 
and Ryan, 2000; Dorrepaal et al., 2009). In contrast, other studies have suggested that warming 
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can lead to a net gain of SOC by promoting the input of biomass into the soil such that it exceeds 
the increase in decomposition (Bond-Lamberty and Thomson, 2010). In most cases, SOC and 
STN increase with decreasing temperature and increasing precipitation (He et al., 2014; England 
et al., 2016; Zhang et al., 2016). SOC and STN also change along elevation gradients, as different 
elevations are characterized by different climatic conditions (Tsui et al., 2013). SOC generally 
increases with increasing elevation (Tsai et al., 2010; Njeru et al., 2017). In addition to climate, 
soil properties, such as pH (Ou et al., 2017), bulk density (Hossain et al., 2015; Chen et al., 2017), 
total phosphorus (Liu et al., 2013; Tian et al., 2017) and soil moisture content (Singh et al., 2011; 
Zhao et al., 2016), are also strongly related to SOC and STN. 

Based on these results, we hypothesize that the SOC and STN values vary among grassland 
types and decrease with decreasing elevation. Furthermore, the concentrations of SOC and STN 
are correlated with meteorological factors and soil properties. Meteorological factors are decisive 
factors in the distributions of SOC and STN, and the values of SOC and STN increase with 
decreasing temperature and increasing precipitation. Our objectives are to (1) estimate the 
concentrations of SOC and STN and their storages and distribution characteristics in different 
grassland types along an elevation gradient, and (2) assess the relationships between temperature, 
precipitation and soil properties and the concentrations of SOC and STN in this arid region. 


2 Materials and methods 


2.1 Site description 


The study area in Fuyun County (45°00'—-48°03'N, 88°10'—90°31'E) is shown in Figures 1 and 2. 
It is located in a border region in the northern part of the Xinjiang Uygur Autonomous Region of 
China and extends from the southern foot of the Altai Mountains to northern margin of the 
Junggar Basin. The study area is 32,186 km?, with an elevation ranging from 317 to 3863 m. It 
lies in the semi-arid and arid climatic belt of the temperate climate zone, and has an annual mean 
temperature of 4.60°C (2007-2016) and a mean annual precipitation of 208.41 mm (2007—2016). 
The Irtysh River and the Ulungur River originate from the Altai Mountains and flow through 
Fuyun County. 
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Fig. 1 Cross sectional diagram of the mountain-basin system landscape 


It encompasses complex landforms alternating between mountains and basins (Figs. 1 and 2). 
Grassland is the main ecosystem type in this area, covering 87.46% of the total area and grazing 
is the primary use of grasslands. There are ten types of grassland in Fuyun County according to 
the grassland resource classification system in China: alpine steppe, alpine meadow, mountain 
meadow, temperate meadow steppe, temperate typical steppe, temperate steppe desert, temperate 
desert steppe, lowland meadow, temperate desert, and swamp. Nomads migrate among different 
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grasslands at different elevations according to the seasons, to utilize grasslands with different 
conditions. 
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Fig. 2 Location of the study area and plots in the mountain-basin system. AM, alpine meadow; MM, mountain 
meadow; TTS, temperate typical steppe; TSD, temperate steppe desert; TDS, temperate desert steppe. 


2.2 Sampling and measurement methods 


We investigated the grassland ecosystem in July 2016 and established a transect from north to 
south. Our study was conducted along this transect across five natural grassland types: alpine 
meadow (AM), mountain meadow (MM), temperate typical steppe (TTS), temperate steppe desert 
(TSD), and temperate desert steppe (TDS) (Fig. 2). The spatial distances among the plots in AM, 
MM, TTS, TSD, and TDS were 12.47, 68.73, 39.31, and 38.61 km, respectively. The plots in AM, 
MM, TTS, TSD, and TDS were distributed along the elevation gradient, which decreased from 
2432 to 772 m. Four sampling plots (10 mx10 m, with an interval of 500—1000 m) were 
investigated in each grassland type. In each plot, three soil profiles with an S-shaped pattern were 
selected, and three quadrats of 1 mx1 m were chosen to investigate the vegetation characteristics. 
We collected core samples from all depth intervals using stainless steel tubes with a volume of 
100 cm? (5.1 cm height and 5.0 cm diameter), to quantify bulk density. Soil samples were 
collected at 0-10, 10-20 and 20-40 cm layers, respectively. We recorded the vegetation species 
and the coverage of each species, and then harvested the aboveground biomass in each quadrat. 
All aboveground parts (green parts and litter) of individual species were cut, collected, and placed 
into envelopes. We measured the biomass aboveground by weighing the aboveground parts of the 
plants after drying at 65°C to a constant weight. The details of the plots are shown in Table 1. 

The temperature and precipitation data were obtained from the China Meteorological Data 
Sharing Service System (http://cdc.nmic.cn/home.do). To overcome the effects of extreme values 
in a single year and reflect the mean temperature and precipitation values from a long time series, 
we used the annual temperature and annual precipitation of 2007—2016 to calculate the ten-year 
mean annual temperature (MAT, in °C) and mean annual precipitation (MAP, in mm) for the 
meteorological stations in Xinjiang Uygur Autonomous Region. To improve the accuracy of the 
calculation results, we referred to the regression tree method (Jing et al., 2013), in which the 
latitudinal and longitudinal coordinates, and topographic elements are added to the climate data 
calculation. Therefore, we used the longitude (Lng, in °), latitude (Lat, in °) and altitude (Alt, in m) 
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Table 1 General conditions of the sampling plots 
Item AM MM TTS TSD TDS 
Lng (°) 89.45 89.41 89.75 89.28 89.00 
Lat (°) 47.61 47.57 46.97 46.84 46.55 
Alt (m) 2432 2066 1414 985 772 
MAP (mm) 394.15 359.42 277.06 227.32 209.80 
MAT (°C) -3.55 -1.71 2.23 5.01 5.88 
Zonal soil Alpine meadow soil Mountain meadow Chestnut soil Brown calcic soil Grey brown 
soil desert soil 
Grazing season Summer Summer Spring/Autumn Winter Winter 
a 65-85 75-90 30-60 10-30 15-30 
Aboveground 129.11+8.21 160.35+7.23 108.56+9.54 76.9+8.65 73.55+10.78 
biomass (g/m*) 
Seriphidium 


Poa alpine 
Festuca ovina subsp. 
Kobresia bellardii 
Carex stenocarpa 
Polygonum viviparum 


Dominant species 


Trifolium 
incarnatum 
Poa angustifolia 
Carex buekii 
Alchemilla pinguis 


Stipa capillata 

Festuca ovina 

Spiraea media 
Caragana pumila 


gracilescens 
Stipa glareosa 
Seriphidium 
terraralbae 
Ceratoides lateens 


Stipa sareptana 
Seriphidium 
borotalense 

Stipa glareosa 

Kochia prostrate 


Note: Lng, longitude; Lat, latitude; Alt, altitude; MAP, mean annual precipitation; MAT, mean annual temperature; AM, alpine meadow; 
MM, mountain meadow; TTS, temperate typical steppe; TSD, temperate steppe desert; TDS, temperate desert steppe. Mean+SE. 


of 54 meteorological stations (with elevations ranging from 34.5 to 3504.4 m) in the Xinjiang 
Uygur Autonomous Region to establish a regression equation with temperature and precipitation. 
We then employed the regression equation to estimate the MAT and MAP for the five plots. The 
MAT and MAP values were calculated using the following formulas: 
MAT= —0.039xLng—1.1095xLat—0.005xAIt+64.204 (R?=0.848), (1) 
MAP= —6.505xLng+31.527xLat+0.093xA It-750.872 (R?=0.78 1). (2) 
SOC and STN were determined using a TOC-Analyzer (Multi N/C 3100, Analytik Jena AG, 
Germany) and the modified Kjeldahl method (Bao, 2005), respectively. The soil available N (AN) 
concentration was measured via the alkali diffusion method (Bao, 2005). Soil total phosphorus 
(STP) was fused with sodium hydroxide, and the concentration was determined via Mo-Sb 
colorimetry. Soil pH (pH) was determined at a soil-water ratio of 1:2.5 (w:v) with a pH meter 
(SG2, Mettler Toledo, Switzerland). Soil bulk density (BD) was measured with a soil wreath knife 
(Bao, 2005). The soil moisture content (SMC) was measured via the oven method. SOC storage 
(t/hm?) and STN storage (t/hm?) were calculated as follows (Xie et al., 2007): 


SOC storage=)"" (1- 0%)SOC, x B, x D ,x0.1, (3) 


STN storage=))" (1—0%)STN, x B, x D;x0.1, (4) 


where SOC;, STN;, Bi and D; are the SOC concentration (g/kg), STN concentration (g/kg), bulk 
density (g/cm), and thickness (cm) of the i layer, respectively. The percentage of rock 
fragments >2 mm was 0%. The total volume of rock fragments accounted for far less than 10% 
for all samples and had a negligible effect in the calculations. 


2.3 Statistical analysis 


The data followed a normal distribution based on the normality test. One-way analysis of variance 
(ANOVA) and the least significant difference (LSD) test were used to analyze the differences in soil 
properties among the five grassland types. The analyses were performed using SPSS, version 20.0 
(IBM, Chicago, IL, USA). The correlations between SOC and STN, meteorological factors and soil 
properties were calculated using Pearson's correlation, and graphs were drawn using R (version 
3.2.2). Redundancy analysis (RDA) is an alternative to canonical correlation analysis and allows the 
relationships between two sets of X and Y variables to be examined. Statistical significance was 
tested via the Monte Carlo permutation method. We used the RDA method to investigate the 
proportion of the variability in SOC and STN that can be explained by environment variables. 
The eigenvalues are proportional to the total variance explained by each axis and which are 
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extracted from soil variables as linear combinations of environment attributes (Lepš and Smilauer, 
1988). RDA was performed using Canoco 4.5 software. In addition, we applied a general linear 
model (GLM) to assess the integrative effects of meteorological factors on the distribution of SOC 
and STN concentrations, which was performed in SPSS, version 20.0 (IBM, Chicago, IL, USA). 


3 Results 


3.1 Concentrations and storages of SOC and STN in different grassland types 


Figure 3 shows the concentrations of SOC and STN in different soil layers in different grassland 
types. The SOC concentrations in each soil layer showed a decreasing trend from AM, MM, TTS, 
TSD to TDS (Fig. 3a). The STN concentrations showed a similar decreasing trend, except for 
MM at the soil depth of 0-10 cm (Fig. 3b). SOC and STN concentrations both decreased with 
depth in all grassland types except TDS (Fig. 3). The SOC and STN concentrations in AM, MM, 
TTS and TSD were significantly higher at the 0—10 cm depth than at the 10—20 and 20-40 cm 
depths (Fig. 3). In the full examined soil layer (0—40 cm), the SOC and STN concentrations 
decreased from AM to TDS among the five grassland types with decreasing elevation, with SOC 
values of 44.77, 35.85, 8.49, 5.89 and 2.26 g/kg (Fig. 4a) and STN values of 4.40, 4.27, 0.99, 0.51 


to 0.15 g/kg, respectively (Fig. 4b). 
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Fig. 3 Soil organic carbon concentrations (SOC, a) and soil total nitrogen concentrations (STN, b) in different 
soil layers in different grassland types. AM, alpine meadow; MM, mountain meadow; TTS, temperate typical 
steppe; TSD, temperate steppe desert; TDS, temperate desert steppe. Bars mean standard errors and the bars 
labeled with different letters (a, b and c for different soil layers in the same grassland type, and A, B, C and D for 
the same soil layer in different grassland types) differ significantly (P<0.05). 
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Fig. 4 Soil organic carbon concentrations (SOC, a) and soil total nitrogen concentrations (STN, b) at a depth of 
0—40 cm in different grassland types. AM, alpine meadow; MM, mountain meadow; TTS, temperate typical 
steppe; TSD, temperate steppe desert; TDS, temperate desert steppe. Bars mean standard errors and the bars 
labeled with different lowercase letters represent significant differences among different grassland types (P<0.05). 


The storages of both SOC and STN decreased from AM to TDS. The SOC and STN storage 
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values were significantly higher in AM and MM than in TTS, TSD and TDS in all soil layers 
(Table 2). In the 0-20 cm layer, the SOC storages of AM, MM, TTS, TSD and TDS accounted for 
56.3%, 60.9%, 60.8%, 59.1% and 43.2%, and STN storages accounted for 58.1%, 63.7%, 59.4%, 
56.5% and 45.9%, respectively, of the total SOC and STN storages in the entire examined soil 
layer (0—40 cm) (Table 2). Therefore, the SOC and STN storage values were greater in the 0—20 
cm layer than in the 20—40 cm layer, except in TDS. 


Table 2 SOC storage and STN storage at different soil depths under different grassland types 


Depth (cm) AM MM TTS TSD TDS 
0-10 65.00+18.07* 43.18+12.05° 14.9643.59° 12.4142.71¢ 3.5441.394 
SOC 10-20 52.66417.97" 46.17412.22" 11.5942.22° 9.9542.72" 3.5841.29° 
(t/hm’) 20-40 91.18430.57" 57.29+10.34° 17.1344.21¢ 15.4846.43¢ 9.35+10.84° 
0-40 208.83443.84* 146.64+27.01° 43.6847.24¢ 37.8446.44° 16.4749.764 
0-10 6.57+1.09° 5.86+1.51° 1.8140.84° 1.0740.19° 0.23+0.06° 
STN 10-20 5.15+1.28° 4.92+1.328 1.26+0.18° 0.8340.19% 0.2740.05° 
(t/hm?) 20-40 8.4443.09 6.1441.51° 2.10+0.74° 1.4640.73¢ 0.59+0.19% 
0-40 20.1644.178 16.9243.16° 5.17+1.65° 3.36+0.79° 1.0940.224 


Note: SOC, soil organic carbon; STN, soil total nitrogen. AM, alpine meadow; MM, mountain meadow; TTS, temperate typical steppe; 
TSD, temperate steppe desert; TDS, temperate desert steppe. The values are the mean+SE and different lowercase letters within the 
same row represent significant differences among different grassland types (P<0.05). 


3.2 Soil properties of different grassland types 


The SMC values of AM, MM, TTS, TSD and TDS were 20.68%, 28.25%, 8.15%, 5.18% and 
2.27%, respectively. The SMC values in the 0-40 cm layer decreased with decreasing elevation 
from AM to TSD, except in MM (Fig. 5a). The BD in the 0—40 cm layer increased from AM to 
TDS with decreasing elevation, except in MM (Fig. 5b). The BD values of AM, MM, TTS, TSD, 
and TDS were 1.25, 1.15, 1.37, 1.70 and 1.74 g/cm’, respectively. The pH values of AM, MM, 
TTS, TSD, and TDS were 4.77, 5.02, 7.94, 8.61, and 8.33, respectively (Fig. 5c). The 
concentration of AN was highest in MM (283.83 mg/kg) and decreased from AM to TTS, TSD 
and TDS (Fig. 5d). The STP concentrations for the five grassland types were 0.48, 0.42, 0.60, 
0.53, and 0.36 g/kg, respectively (Fig. 5e). 


357 (@) a 2.0r (b) 10; (c) ‘ 
30L 8} aa 8 c b 
~ 25| Tle f 
$7 b 514p 6 
S 2f 1.2} zf e ¢ 
S 1.04 rai 
Z 15} Bost 4t 
10| c j 6} > 
Hie i 
0 1 an 0 , , 


AM MM TTS TSD TDS AM MM TTS TSD TDS AM MM TIS TSD TDS 


30 (d) a eo 
“ab D 

= 250 b z 06 

J 200 3 

E 150 5 

z ; 

2 100 d 5 

3 50 d ps 

2° dg Be 

AM MM TTS TSD TDS AM MM TTS TSD TDS 


Fig. 5 Soil properties of different grassland types. Bars mean standard errors and the bars labeled with different 
lowercase letters differ significantly (P<0.05). SMC, soil moisture content; BD, bulk density; pH, pH value; AN, 
available nitrogen; STP, soil total phosphorus; AM, alpine meadow; MM, mountain meadow; TTS, temperate 
typical steppe; TSD, temperate steppe desert; TDS, temperate desert steppe. 
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3.3 Correlations among meteorological factors, soil properties and the concentrations of 
SOC and STN 


Pearson's correlation showed close relationships among the SOC and STN concentrations, 
meteorological factors and soil properties (Fig. 6). The SOC and STN concentrations exhibited 
significant positive correlations with MAP (7=0.933, P<0.01; r=0.938, P<0.01, respectively), AN 
(r=0.905, P<0.01; r=0.948, P<0.01, respectively) and SMC (r=0.875, P<0.01; r=0.922, P<0.01, 
respectively) and negative correlations with MAT (r= —0.928, P<0.01; r= —0.933, P<0.01, 
respectively), pH (r= —0.931, P<0.01; r= —0.946, P<0.01, respectively) and BD (r= —0.809, 
P<0.01; r= —0.846, P<0.01, respectively). The SOC concentrations were also highly correlated 
with the STN concentrations (7=0.977, P<0.01) (Fig. 6). 
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Fig. 6 Correlation diagram of SOC and STN concentrations, meteorological factors, and soil properties 
measured in the 0—40 cm soil layer. The figure shows the correlation coefficients, and the colored solid circles 
represent significant correlations (P<0.05). SOC, soil organic carbon; STN, soil total nitrogen; STP, soil total 
phosphorus; AN, available nitrogen; pH, pH value; BD, bulk density; SMC, soil moisture content; MAT, mean 
annual temperature, MAP, mean annual precipitation. 


Furthermore, the RDA method was employed to assess the factors that best explained the 
variability in the SOC and STN concentrations. Seven environmental factors were selected for 
further RDA analysis after a unified treatment of different types of data. The MAP and MAT 
represent the meteorological factors, and five other factors (AN, STP, BD, pH, and SMC) 
represent the soil properties. The RDA ordination accurately reflected the relationships among 
SOC, STN and the environmental factors, which yielded the canonical coefficients of the 
environmental factors for the two axes (Fig. 7; Table 3). The Monte Carlo permutation test 
showed that the eigenvalues for the first and second axes accounted for 92.5% of the variation in 
SOC and STN (P<0.01; Table 3). Axis 1 was positively correlated with AN, SMC, and MAP and 
negatively correlated with pH, BD, MAT and STP. Axis 2 was positively correlated with AN, 
SMC, MAP and STP, but negatively correlated with the other impact factors (Table 3; Fig. 7). The 
correlations and eigenvalues resulting from the RDA ordination, indicated that the MAP and MAT 
appeared to be the main regulators of the variation in the SOC and STN concentrations (MAP: 
Taxisı1=0.94, and MAT: raxisı= —0.93, P=0.002; Table 3). However, the effects of these two factors 
on the SOC and STN concentrations need to be further discussed. 

Moreover, we applied a GLM to assess the integrative effects of meteorological factors on the 
variation in the SOC and STN concentrations and determined the amount of the variability that 
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could be explained by MAT and MAP separately. The results showed that MAP and MAT together 
explained 97.85% and 98.38% of the overall variation in SOC and STN concentrations, 
respectively, in the 0—40 cm soil layer. MAP, as the most important meteorological factor, 
explained 54.83% and 54.79% of the variation in SOC and STN concentrations, respectively, in 
the 0—40 cm soil layer (Table 4). 
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Fig. 7 Redundancy analysis (RDA) ordination diagram for soil organic carbon (SOC) and soil total nitrogen 
(STN) concentrations and environmental factors for the 0—40 cm soil layer. The numbers 1, 2, 3, 4 represent the 
plots in AM; 5, 6, 7, 8 represent the plots in MM; 9, 10, 11, 12 represent the plots in TTD; 13, 14, 15, 16 represent 
the plots in TSD; and 17, 18, 19, 20 represent the plots in TDS. AN, available nitrogen; pH, pH value; BD, bulk 
density; SMC, soil moisture content; MAT, mean annual temperature; MAP, mean annual precipitation; STP, soil 
total phosphorus. 


Table 3 Canonical coefficients of environmental factors with the first two axes of the redundancy analysis 


Impact factors Axis 1 Axis 2 
AN 0.91 0.18 

pH —0.93 —0.06 

BD -0.81 —0.16 
SMC 0.88 0.21 
MAP 0.94 0.10 
MAT —0.93 —0.01 
STP —0.12 0.02 
Eigenvalues 0.90 0.024 


Cumulative percentage 


h 90.1 92.5 
of variance 


Note: AN, available nitrogen concentration; pH, pH value; BD, bulk density; SMC, soil moisture content; MAT, annual mean 
temperature; MAP, mean annual precipitation; STP, soil total phosphorus. 


Table 4 Effects of MAP and MAT on SOC and STN concentrations in the 0—40 cm layer 


SOC STN 
Parameter 
df MS SS% df MS SS% 
MAP 1 850.56 54.83 1 10.98 54.79 
MAT 1 667.34 43.02 1 8.73 43.59 
Residuals 57 33.42 2.15 57 0.36 1.62 


Note: MAP, mean annual precipitation; MAT, mean annual temperature; SOC, soil organic carbon; STN, soil total nitrogen. df, degrees 
of freedom; MS, mean squares; SS%, proportion of variance explained by the variable; the presented results were obtained from GLM 
analysis. 
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4 Discussion 


4.1 SOC and STN variations among grassland types and soil depths 


In our study, we investigated the concentrations and storages of SOC and STN in different 
grassland types in the MBS of Fuyun County, Xinjiang, China. Estimates of SOC and STN 
storages at the regional scales commonly use vegetation type as a reference (Yimer et al., 2006). 
Vegetation is the primary source of SOC and STN and can significantly affect the SOC and STN 
concentrations through variations in the quantity and quality of the organic matter input into the 
soil, resulting in variations in SOC and STN concentration and storage values (Podwojewski et al., 
2011; Tsui et al., 2013). The grasslands of AM and MM primarily consist of gramineous species 
and leguminous species grasses with relatively high aboveground biomass, representing relatively 
high grassland productivity. In contrast, the grassland of TTS primarily consists of short grasses 
and semi-shrubland species grasses, which feature moderate aboveground biomasses and the 
vegetation of the TSD and TDS grasslands, are primarily semi-shrubland species grasses and 
bunchgrasses with low aboveground biomasses (Table 1). The shifts in aboveground biomass and 
constitutive species from AM to TDS (i.e., from higher elevation to lower elevation) was likely 
responsible for the variation in the C and N inputs from plant biomass to the soils (Wynn et al., 
2006). Therefore, the SOC and STN concentrations and storages decreased from AM to TDS in 
each of the investigated soil layers, except for the STN concentration in the 0—10 cm layer. The 
STN concentration at 0-10 cm layer in MM was higher than that in AM, possibly due to the 
vegetation species composition. In addition to the gramineous species, leguminous species 
(Trifolium incarnatum), which are able to fix N, were present in MM. The SOC and STN storage 
values for 0—40 cm decreased from AM (208.83 and 20.16 t/hm?, respectively), MM (146.64 and 
16.92 t/hm?, respectively), TTS (43.68 and 5.17 t/hm’, respectively), TSD (37.84 and 3.36 t/hm?, 
respectively) to TDS (16.47 and 1.09 t/hm?, respectively) in our study, which is in agreement with 
the overall decreased trend in China (in which the median soil C storage values in alpine meadow, 
alpine steppe, temperate typical steppe, temperate deserted steppe and temperate desert are 182, 
170, 123, 87 and 62 t/hm’, respectively, at soil depths of 0-100 cm (Ni, 2001). Yang et al. (2007) 
estimated the SOC storages in five biomes in China and showed the same decreased trend 
(meadow>steppe>desert). Liu et al. (2012) observed that the SOC and STN storage values in 
alpine meadow (87.0 and 8.1 t/hm?, respectively) were higher than in desert steppe (70.9 and 9.9 
t/hm?, respectively) and desert (43.9 and 6.8 t/hm?, respectively) in the upper 1-m soil layer on the 
QTP. However, due to spatial heterogeneity in the soil, the SOC and STN storage values in our 
studied grassland types were different than those in other regions. Although the grassland types 
were between our study and other studies, the soil types and climate conditions were very 
different, both of which play important roles in determining the SOC storage distributions (Liu et 
al., 2012). 

Our results showed that, except in TDS, the SOC and STN concentrations decreased with depth 
(Fig. 3), and that total SOC and STN storage values in the 0—10 and 10-20 cm soil layers were 
higher than in the 20—40 cm layer (Table 2). This finding is in agreement with other observations 
made on the QTP (Yang et al., 2010) and in other ecosystems (Zhang et al., 2013; Dorji et al., 
2014). Topsoil is more readily affected by litter falling on the soil surface. The C and N inputs to 
the subsoil primarily depend on the migration of C and N from the topsoil to subsoil. In our study 
area, AM and MM usually exhibited thick humus, which is considered to represent the most 
important soil C sink (Velichko et al., 2010). However, this situation was not observed in TDS. 
The 20—40 cm soil layer in TDS showed the higher SOC and STN concentrations and storages 
values. Liu et al. (2012) previously observed the same trend in the desert steppe of the QTP. 
These results might be related to the scarcity of plants and low plant coverage on the surface in 
TDS, resulting in a low C input. Additionally, SOM decomposition in the topsoil is enhanced by 
the relatively higher temperatures in the topsoil than in the subsoil. The low plant coverage and 
high temperature also result in less soil moisture in the topsoil. Some studies have shown that low 
soil moisture was advantageous for the allocation of biomass belowground (Gansert, 1995). 
Another possible reason is the frequently observed high sediment input, where fine sand and silts 
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have been deposited on surface, frequently burying the humic horizons of earlier soil formation 
beneath fresh sediment (Liu et al., 2012). 


4.2 SOC and STN relationships with meteorological factors 


The concentrations and storages of SOC and STN decreased with decreasing elevation throughout 
the entire 40 cm soil profile (Fig. 4; Table 2), similar to the findings of previous studies in other 
regions (Garten and Hanson, 2006; Podwojewski et al., 2011). Previous study showed that MAP 
was positively correlated with SOC and STN concentrations, while MAT was negatively 
correlated. Our results are also consistent with those of previous studies conducted at a global 
scale (Ganuza and Almendros, 2003; Wang et al., 2010). However, there is a lack of consensus at 
the regional scale. Yang et al. (2010) found that SOC concentration increased with MAT in 
Tibetan grasslands, while our study showed that the SOC and STN concentrations decreased with 
MAT, but were significantly and positively correlated with MAP (Figs. 6 and 7). In Fuyun County, 
differences in temperature and precipitation were observed among different elevations. Our 
results indicated that the variations in MAP and MAT along the elevation gradient controlled the 
SOC and STN concentrations (Fig. 7; Table 3). 

The temperature and precipitation affect the SOM input from different vegetation species litter 
production, as well as the output through decomposition and mineralization (Giardina and Ryan, 
2000). At high elevations, higher levels of precipitation can result in an increase in plant biomass, 
thus further resulting in greater organic matter and nitrogen input. In addition, the decomposition 
and mineralization of SOM are inhibited in low temperature environments, resulting in the 
accumulation of soil C and N (Tsai et al., 2010). In contrast, in low-elevation areas in arid regions, 
the higher temperature reduces the water use efficiency of plants by increasing evapotranspiration 
(Lal, 2004). Hence, the water deficit induced by higher temperatures and lower precipitation may 
restrict plant growth. Furthermore, the higher decomposition and mineralization rates resulting 
from the higher temperatures may accelerate the decrease in SOC and STN concentrations 
(Simmons et al., 1996). This effect was significant in Fuyun County, an arid region, since 
vegetation growth is mainly limited by natural precipitation (Liu et al., 2013). The GLM analysis 
suggested that precipitation explained a higher proportion of the distribution of SOC and STN 
concentrations than temperature (Table 4). These results imply that precipitation is a limiting 
factor for vegetation growth and vegetation productivity in the examined arid grassland 
ecosystem, as a small increase in the precipitation may significantly increase the productivity of 
vegetation, thus contributing to the accumulation of SOC and STN concentrations (Jobbágy and 
Jackson, 2000; Callesen et al., 2003). The results showed that the differences in the climate 
conditions, such as temperature and precipitation, caused by the differences in elevation are the 
most important environmental factors affecting the spatial heterogeneity of SOC and STN. 

To simply estimate SOC and STN concentrations in arid regions with a MBS structure, we 
established predictive models for SOC and STN concentrations based on MAP (in mm) and MAT 
(in °C). The predictive models for the SOC and STN concentrations are as follows: 

SOC=35.685xMAT+2.054xMAP-639.562, (5) 
STN=4.081xMAT+0.233xMAP-—72.839. (6) 

The equations effectively fit the relationships between the meteorological factors and the SOC 

and STN concentrations, with R? values of 0.904 and 0.919, respectively. 


4.3 SOC and STN relationships with soil properties 


The C cycle is closely linked with the N cycle through production and decomposition (Liu, 2001). 
Our results also showed that the STN concentration presented a significant positive correlation 
with the SOC concentration in the 0—40 cm soil layer of the different grassland types (Fig. 6) 
because the main N sources were litter, animal residue, biological nitrogen fixation, and organic 
matter (Bai et al., 2001). 

AN was found to be positively correlated with SOC and STN in the present study (Fig. 6). The 
AN included the inorganic N derived via the mineralization of SOM (Hodge et al., 2000) and 
organic N derived via the decomposition of SOM (Kielland, 1994). Therefore, a high SOM 
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concentration usually is associated with high SOC, STN, and AN concentrations. 

Our results indicated that SMC was positively correlated with SOC and STN (Fig. 6). This 
conclusion agrees well with the findings of previous studies (Sun and Wang, 2016; Wang et al., 
2016). The anoxic decomposition of SOM tends to be inhibited under higher SMC conditions, 
resulting in the accumulation of SOC. Moreover, SMC also affects N, as a higher SMC can limit 
soil microbial activities, producing an environment unfavorable to the mineralization and 
decomposition of soil organic nitrogen (Bai et al., 2004; Zhao et al., 2016). Therefore, a high 
SMC can lead to a high STN concentration. In our study, the SMC in the 0—40 cm layer decreased 
from AM to TDS with decreasing elevation, with the exception of MM (Fig. 5a). MM might not 
follow the general trend because it was located in a valley and the Irtysh River flowed through the 
sites. The slope of the riverbed was small, the riverbed was broad, the runoff was large, and the 
current velocity was low. In addition, the underground water level was high. Therefore, the SMC 
was higher in MM. 

We found that soil pH exhibited a significant negative correlation with SOC and STN (Fig. 6). 
This finding is in agreement with the conclusions of other studies (Njeru et al., 2017; Ou et al., 
2017). Soil pH influences SOC and STN by regulating the microbial structure, function and 
activities (Bai et al., 2002). The decomposition of SOM is inhibited at lower pH values and 
accelerated at higher pH values (Motavalli et al., 1995). The SOC, STN, and pH interact with 
each other via meteorological factors. Cui et al. (2005) reported that topsoil pH is closely 
correlated with precipitation at a large spatial scale in arid areas. Our results indicated that the pH 
was significantly negatively correlated with MAP (Fig. 6). Hence, high precipitation resulted in 
increases in SOC and STN (Figs. 6 and 7), as well as in the root systems of plants. Consequently, 
more carbonic and organic acids were produced by root exudation and soil respiration, which 
could be responsible for the decrease in pH. 

BD was significantly negatively correlated with SOC and STN concentrations in the 0—40 cm 
soil layer (Fig. 6). This relationship has also been recognized by other studies (Zhang et al., 2016; 
Njeru et al., 2017). The SOC and STN concentrations were negatively correlated with BD, 
possibly because the mineralization of SOC and the nitrification of ammonium nitrogen were 
suppressed in soils with high BD values (De Neve and Hofman, 2000). A low BD soil can store 
more SOC and STN, as SOC and STN can be mobilized in porous spaces within the soil matrix 
(Anh et al., 2014). Furthermore, higher levels of biomass can produce more soil nutrients via 
decomposition processes, and the soil BD was negatively correlated with soil nutrients. Lister et 
al. (2004) also reported that BD can be influenced by root development and soil microbial 
activities and would decrease under higher nutrient levels because of the creation of more and 
larger soil pores. In addition, the BD in the 0-40 cm layer increased from AM to TDS with 
decreasing elevation, except in MM (Fig. 5b). The aboveground biomass was highest in MM, and 
there was a thick layer of humus in the surface soil due to the abundance of litter and roots, which 
resulted in a loose and porous structure (and thus a lower BD) for the soil in MM. 

STP was not correlated with SOC and STN in our study (Fig. 6), similar to the findings of 
previous studies in other regions (Gao et al., 2016; Hu et al., 2018). P is primarily derived from 
rocks and its poor mobility results in a relatively high independence from other nutrients. Hence, 
its spatial distribution is highly heterogeneous. The factors that influence the STP patterns are 
usually quite complex. In addition to the soil parent material, the climate and biota also play 
important roles (Gao et al., 2016). Although STP concentrations have been reported to associated 
with SOM (Vincent et al., 2010; Camargo et al., 2013), the STP concentrations in our study 
showed no significant correlation with the SOC or STN concentrations. 

Our results accurately reflect the characteristics of the concentration and storage values of SOC 
and STN in different grassland types along an elevation gradient, as well as their relationships 
with meteorological factors and soil properties. The present study is an effective complement to 
the available data on regional SOC and STN concentrations and storages estimates in grasslands 
in arid regions with a MBS structure. 
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5 Conclusions 


In Fuyun County, SOC and STN were found to vary significantly with the grassland type. Our 
results indicated that the SOC and STN concentration and storage values in the 0—40 cm soil layer 
decreased in the five investigated grassland types along a decreasing elevation gradient and also 
decreased with increasing soil depth, except in TDS. The SOC and STN concentration and storage 
values of different grassland types in Fuyun County presented obvious vertical zonality 
characteristics. The results of Pearson's correlation and RDA analysis showed that the SOC and 
STN concentrations in the 0—40 cm soil layer were significantly positively correlated with MAP, 
SMC and AN and negatively correlated with MAT, pH and BD. The SOC and STN concentrations 
were significantly affected by meteorological factors. MAP and MAT together explained 97.85% 
and 98.38% of the overall variations in the SOC and STN concentrations, respectively, in the 
0—40 cm soil layer. Moreover, MAP explained 54.83% and 54.79% of the variations in the SOC 
and STN concentrations in the 0—40 cm soil depth, respectively. Our results demonstrated that 
precipitation and temperature are simple and effective predictors of the SOC and STN 
concentration and storage values in grasslands in arid regions with a MBS structure. 
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